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Abstract. In view of the fact that the data on neutrino mixing are still compatible with a 
situation where Bimaximal mixing is valid in first approximation and it is then corrected by 
terms of order of the Cabibbo angle, we present examples where these properties are naturally 
realized. The models are supersymmetric in 4-dimensions and based on the discrete non-Abelian 
flavour symmetry S4. 



1. The flavour puzzle and the neutrino mixing patterns 

The information on masses and mixings are encoded in the Yukawa terms, which are theoretically 
unpredicted. On the contrary, from the experiments we know that the charged fermions have 
strong hierarchical masses, spanning from tenths of MeV up to few hundreds GeV, while the 
neutrinos do not show a strong hierarchy and only upper bounds of few eV have been put on 
their masses. Furthermore, the solar and the atmospheric anomalies find an elegant explanation 
in the oscillation of three active massive neutrinos and the corresponding frequencies have been 
measured in terms of mass squared difference of neutrinos eigenvalues. In table [TJ we can see 
the result of two independent global fits on the neutrino oscillation data. An uncertainty is still 
present on the sign of the atmospheric mass squared difference and this reflects our ignorance 
on the type of the neutrino spectrum: normal hierarchy, inverse hierarchy or quasi degeneracy. 

When looking at the mixings, other differences are underlined: the quark mixing matrix, 
Vckm, has only small angles which can be expressed in terms of powers of the Cabibbo angle, 
A ~ 0.23, through the Wolfenstein parametrisation, while the lepton mixing matrix, Upmns, 
presents two very large angles and the third compatible to be vanishing: in particular the 
atmospheric angle, in the 23 sector, is compatible with the maximal value well inside the la 
region; the solar angle, in the 12 sector, is large, but about 5cr's from being maximal; the reactor 
angle, in the 13 sector, is compatible to be vanishing at about la level as a result of global fits 
on the neutrino oscillation data, but only an upper bound of about A has been put on sin #13 
by reactor experiments. A very good approximation to the lepton mixing matrix is provided by 
the so-called Tri-Bimaximal (TB) patter [3111], 



sin 2 6»i2 = 1/3 



sin 2 23 = 1/2 



sin 6*i3 = 




Table 1. Neutrino oscillation parameters from independent global fits [HE]- 
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which agrees at the la level with the data. Note that U does not depend on the mass 
eigenvalues, in contrast with the quark sector, where the entries of the CKM matrix can be 
written in terms of the ratio of the quark masses. Moreover it is a completely real matrix, since 
the factors with the Dirac phase vanish (the Majorana phases can be factorized outside). The 
best measured neutrino mixing angle 9\2 is just about la below the TB value, while the other 
two angles are well inside the la interval. 

In a series of papers [51 El [H |9] it has been pointed out that a broken flavour symmetry 
based on the discrete group A4 appears to be particularly suitable to reproduce this specific 
mixing pattern as a first approximation. Other solutions based on alternative discrete or contin- 
uous flavour groups have also been considered [TOl [Til 021 US E] , but the A± models have a very 
economical and attractive structure, e.g. in terms of group representations and of field content. 
In all these models, when the symmetry is broken, some corrections to the mixing angles are 
introduced: in general all of them are of the order of Xq and therefore these models indicate 
a value for the reactor angle which is well compatible with zero (for a different approach see |15j). 

There is an experimental hint for a non- vanishing reactor angle [H[2] and, if a value close to the 
present upper bound is found in the future experiments, this could be interpreted as an indication 
that the agreement with the TB mixing is only accidental. Looking for an alternative leading 
principle, it is interesting to note that the data suggest a numerical relationship between the 
lepton and the quark sectors, known as the complementarity relation, for which #12 + Ac — vr/4 
|16| [TT] 118] . However, there is no compelling model which manages to get this nice feature, 

without parameter fixing. Our proposal is to relax this relationship. Noting that yj m^/rrir ~ Ac, 
we can write the following expression, which we call weak complementarity relation 




The idea is first to get a maximal value both for the solar and the atmospheric angles and then 
to correct #12 with relatively large terms. To reach this task, the bimaximal (BM) pattern |19j 
can be extremely useful: it corresponds to the requirement that #13 = and #23 = #12 = tt/4. 
The unitary matrix which corresponds to this patter is given by 

/ 1/V2 -1/V2 \ 

jjBM = y 2 x j 2 -l/y/2 . (3) 

\ 1/2 1/2 +1/V2 j 



Note that U BM does not depend on the mass eigenvalues and is completely real, like the TB 
pattern. The most general mass matrix, m BM , diagonalized by this mixing is \i — r symmetric 
and satisfies to an additional symmetry for which (m^ )i 1 = {m„ M )2,2 + ( m ^ M )2,3'- 



m^M = \ y z x - z \ (4) 




Starting from the BM scheme, the corrections introduced from the symmetry breaking must 
have a precise pattern: 5 sin 2 #12 ^ Ac, while <5sin 2 #23 < an d <5sin#i3 < Ac in order to be 
in agreement with the experimental data. This feature is not trivially achievable. 



2. The model building 

In this part we present two flavour models \20\ [2T] in which the neutrino mixing matrix at the 
leading order (LO) is the BM scheme, while the charged lepton mass matrix is diagonal with 
hierarchical entries; moreover the models allow for corrections that bring the mixing angles in 
agreement with the data. The strategy is to use the flavour group Gf = S4 x Z4 x U(1)fn, 
where 54 is the group of the permutations of four objects, and to let the SM fields transform 
non-trivially under Gf; moreover some new fields, the flavons, are introduced which are scalars 
under the SM gauge symmetry, but transform under Gf] these flavons, getting non-vanishing 
vacuum expectation values (VEVs), spontaneously break the symmetry in such a way that two 
subgroups are preserved, Gg = Z4 in the charged lepton sector and G v = Z2 x Z^ in the neutrino 
sector. It is this breaking chain of S4 which assures that the LO neutrino mixing matrix, U v , 
is the BM pattern in the basis of diagonal charged lepton mass matrix. The additional terms 
Z4 x U(1)fn forbid dangerous operators and allow for the correct charged lepton mass hierarchy. 

The first model |20| deals only with the lepton sector in a supersymmetric SM scenario, 
implementing the weak complementarity. In the second model |21j . we extend also to the quark 
sector in a Pati-Salam GUT context. 



2.1. The lepton model 

We formulate our model in the framework of the See-Saw mechanism. For this we choose the 3 
generations of left-handed lepton doublets I and of RH neutrinos v c transforming as (3, 1) under 
S4 x while the RH charged leptons e c , \i c and r c transform as (1, —1), (1', —%) and (1, —i), 
respectively. The S4 symmetry is then broken by suitable flavons, ipt and Xli transforming as 
(3,i) and (3',i), and ip„ and £„, transforming as (3,1) and (1,1). A flavon 9, carrying only a 
negative unit of the U(1)fn charge, acquires a VEV and breaks U(1)fn- In view of a possible 
GUT extension of the model at a later stage, we adopt a supersymmetric context, so that two 
Higgs doublets H u ^, invariant under 54 x Z4, are present in the model. The usual continuous 
U(1)r symmetry, related to i?-parity, is implemented in the model. Supersymmetry also helps 
producing and maintaining the hierarchy {H u ^) = v u d <C Ay where Af is the cutoff scale of the 
theory. 

The complete superpotential can be written as 

W = W e + W u + W d , (5) 

where w^ is responsible for the flavon VEV alignment as discussed in |20j , while w e and w v refer 
to the charged lepton and neutrino sectors and can be written as 

we = y ^r^e c {^m)H d + ^^e%l X tXt)H d + ^^e c {l mX t)H d + 

+ ^jj^XiyH d + ^T c (iip e )H d + ... ' ' (6) 
w v = y(v c £)H u + Mk s {v c v c ) + a(v c v c t u ) + b{y c v c ^ v ) + ... 



indicating with (. . .) the singlet 1, with (. . .)' the singlet 1' and with (. . .)r the representation 
R (R = 2, 3, 3'). Note that the parameter M defined above is dimensionless. In the above 
expression for the superpotential w, only the lowest order operators in an expansion in powers of 
1/Aj are explicitly shown. Dots stand for higher dimensional operators. The stated symmetries 
ensure that, for the leading terms, the flavons that appear in w e cannot contribute to w u and 
viceversa. 

We showed in |20] that the potential corresponding to wa possesses an isolated minimum for 
the following VEV configuration: 






(<p t )/A f ~\ 1 \ v, (xi)/A f ~\ \cv, (^)/A/~ 1 \v', (U/A/~cV, (7) 



where c and d are order one coefficients which parametrize the ratio among the different VEVs. 
Similarly, the Froggatt-Nielsen flavon 9 gets a VEV, determined by the D-term associated to 
the local U(1)fn symmetry, and it is denoted by (9}/Af = t. 

With this VEVs configuration, the charged lepton mass matrix is diagonal 

M e ~diag(t 2 , t, l)vv d /y/2, (8) 

where for simplicity we omit to indicate all the coupling constants and the coefficient c. As a 
result, at the LO, there is no contribution to the lepton mixing matrix from the diagonalisation 
of charged lepton masses. It is possible to estimate the values of v and t by looking at the mass 
ratios of charged leptons: 

m M /m r ~i, m e /m^~ v t. (9) 

In order to fit these relations with the data, we must have approximately t ~ 0.06 and v ~ 0.08 
(modulo coefficients of 0(1)). 

The light neutrino Majorana mass matrix, given by the See-Saw relation, is diagonalised by 
the BM matrix and the eigenvalues are 

I y^ 1 I li^ 

U^m v U u = diag(mi, m 2 , m 2 ) with m, = , (10) 

where Mi are the masses of the heavy RH neutrinos 

Mi = 2 M + v'(a- V26c') A/ , M 2 = 2 M + v' (a + V2bc') A/ , M 3 = 2 \M + av'\ A f . 

(11) 

Notice that since the neutrino sector is not charged under the symmetry, we have operators 
of dimension 5 which contribute to the neutrino masses and may correspond to some heavy 
exchange other than the right-handed neutrinos v c . When considering the interesting domain 
of parameters, we find that this effective contribution is subdominant. 

At the LO, the light neutrino mass matrix depends on only 2 effective parameters indeed the 
terms M and av' enter the mass matrix in the combination F = M + av' . The coefficients y£\ 
Ufj,, Ut, V, o, and b are all expected to be of 0(1). A priori M could be of 0(1), corresponding 
to a RH neutrino Majorana mass of 0(A/), but, actually, we saw that it must be of the same 
order as v' . 

To summarise, at the LO we have diagonal and hierarchical charged leptons together with 
the exact BM mixing for neutrinos. It is clear that substantial the next-to-leading order (NLO) 
corrections are needed to bring the model to agree with the data on Oyi- A crucial feature of 



our model is that the neutrino sector flavons (p„ and are invariant under Z4 which is not 
the case for the charged lepton sector flavons (p$ and xe- The consequence is that (p u and £„ 
can contribute at the NLO to the corrections in the charged lepton sector, while at the NLO 
<pi and xt cannot modify the neutrino sector couplings. As a results the dominant genuine 
corrections to the BM mixing only occur at the NLO through the diagonalisation of the charged 
leptons. Without entering in the details of the NLO discussion (see [20]); we find that the NLO 
corrections, coming from the higher order terms, are not democratic and the final corrected 
mixing angles are 

sin 2 6>i2 = - - _/|( c l +C2W , sin 2 #23 = ~ , sin6»i 3 = _( Cl - c 2 )v' , (12) 

where c\ and C2 are coefficients of order one. When v' is of the order of the Cabibbo angle, then 
#12 is brought in agreement with the experimental data; in the meantime the reactor angle is 
corrected of the same amount, suggesting a value for #13 close to its present upper bound. Note 
that the atmospheric angle remains uncorrected at this order. Any quantitative estimates are 
clearly affected by large uncertainties due to the presence of unknown parameters of order one, 
as we can see in figure [TJ but in our model a value of #13 much smaller than the present upper 
bound would be unnatural. 
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Figure 1 . sin #13 as a 

function of sin 2 612 is plotted, 
following eqs. (12). The 
parameters c\ and C2 are 
treated as random complex 
numbers of absolute value 
between and 0.30. The gray 
bands represents the regions 
excluded by the experimental 
data [IJ[2]: the horizontal one 
corresponds to the 3<r-upper 
bound for sin 2 #13 of 0.46 and 
the vertical ones to the region 
outside the 3<r error range 
[0.26- 0.37] for sin 2 6»i 2 . 



It is then interesting to verify the agreement of the model with other sectors of the neutrino 
physics, such as the 0^2/3-decay and the leptogenesis (See [22j [23] for a general approach). The 
result of the analysis is that the model presents a normal ordered - moderate hierarchical or quasi 
degenerate - spectrum with a suggested lower bound for the lightest neutrino mass and for the 
effective 0^2/3-mass parameter \m ee \ of about 0.1 meV. On the other hand it is compatible with 
the constraints from leptogenesis as an explanation of the baryon asymmetry in the Universe. 

2.2. The full model 

In this section we are interested in the extension to the quark sector. A first attempt is to 
adopt for quarks the same representations under S4 that have been used for leptons: the left- 
handed quark doublets q transform as a triplet 3, while the right-handed quarks (ii c , d c ), (c c , s c ) 
and (t c , b c ) transform as 1, 1' and 1, respectively. We can similarly extend to quarks the 
transformations of Z\ (and U(1)r) given for leptons. As a result, it is easy to see that the 
quark mass matrices are diagonal, at the LO on the expansion parameters, exactly as for the 



charged leptons and to account for the correct mass hierarchies the U(1)fn has to be suitably 
implemented. At this level the CKM matrix is the unity matrix and to get realistic mixings, 
the higher-order corrections should switch on off-diagonal entries with a well-defined pattern: 
(12) ~ A, (23) ~ A 2 and (13) ~ A 3 . By an explicit computation we find the following result for 
the quark mass matrices: in terms of order of magnitude 

v v' t 2 \ / v t 3 vv' t 3 v v' t 3 

,,/2 + 1 vv d )l/T _ I +2 4.2 „J2+2 ^ VV ' 1 ' 





M d = \ v't t v' 2 t M u = \ v't 2 t 2 v' 2 t 2 \ -^. (13) 




v 12 



V2 • 



Calculating now the unitary matrices which diagonalise M* d M d and MlM u we find 

/ 1 v' v' \ 
V d ~V u =\ -v> 1 v' 2 . (14) 
V -v' -v' 2 1 / 

At a first sight, barring cancellations among the single entries, we can see that the CKM matrix 
V = V^Vd should be similar to V u and and as a consequence it cannot correctly describe the 
quark mixings: while the entries (12) and (23) well reproduce the measured values, the large 
values in the (13) entry would require a large fine-tuning of order A 2 . 

An alternative possibility is to further investigate on the complementarity relations: 



'12 



+ A-7T/4, 023 + A 2 ~ -tt/4 . (15) 



These equations suggest that the angles in the CKM and PMNS matrices may have a common 
origin which can be motivated for example in Pati-Salam models, where the following relation 
holds, 

U e ~ V d . (16) 

We can use this to write the CKM and PMNS matrices as 

U = R23 (~f ) -Rl3(A)i?12 (f - A) = (^23(1)^13(^12(^)^12(1) 



Ue U v 



V = R 12 (X) = (i? 2 3(^) J Ri3(A)^i 2 (A)) t J R 23 (^)i?i3(A)i?i2(A)). 



(17) 



v u v d 



Here Rij{a) stand for rotations in the (ij) plane of the angle a (apart from coefficients 0(1) 
in front of each angles). The coefficients of the angles in the rotations in and V d should be 
such that the rotations cancel each other in the (13) sector, but not in the (12) sector. Thus 
we should introduce terms which distinguish between the up- and down-quark sectors, indeed 
possible within the Pati-Salam context. 

Moving to the explicit form of the mass matrices, the generic Majorana neutrino mass matrix 
m u which is diagonalised by U u = R\2 (f) , 

m f-9 = Rl2 ^y muRl2 ^j t ( 1 8) 

is given by 

m v ~ I b a I . (19) 




Considering the charged lepton mass matrix M e , the product Aft M e should be diagonalised by 
the action of Vd as in eq. 17 



^i2(-A) J R 13 (-A)i? 23 f-*J Aft M e i? 23 R 13 {X)R 12 (X) . (20) 
In the limit m e — )• we find the generic structure for the product Aft M e : 



AftM e ~ *f (A 1 1 J +^ I A^ 1 -1 ] + 0(A 2 ) , (21) 



A A \ _2/0 A -A 



that can be obtained if Af e is given by 



\ / 

M e ~^| +^ A 1 -1 |+0(A 2 ). (22) 
v 2 \ A 1 1 / v 2 \ o o o 



It is interesting to note that, moving to the basis of diagonal charged leptons and considering 
only the LO terms, the neutrino mass matrix results to be of the classical BM type. 

the relation M e ~ Af^ follows and therefore the down-quark matrix has a 
Looking at eq. (17) we find that also Af u should have a simi- 



From eq. 16 



similar structure as in eq. 22 



lar structure. We find that we can satisfy the constraint on the eqs. (12) and (13) rotations if 
the third columns of M u and are proportional to each other, but the second columns are not. 



So far we have not given yet any explanation of the origin of these mass matrices and mixings 
and we leave this analysis to [21J. Here we only say that such a construction is possible in a 
Pati-Salam realisation where the flavour symmetry is S4 x Z4 x U(1)fn- as f° r the non-GUT 
model described above, key points are a suitable choice of the group representations for the 
flavons and the particular VEV misalignment whose effects are a reactor angle and a deviation 
from 7r/4 of the solar angle of the order of A, while introducing small deviations of the order A 2 
from the maximal value of the atmospheric angle. 

A difficulty with respect the non-GUT model refers to the study of the gauge coupling running 
and of the Higgs potential: while in a general Pati-Salam model, in particular without any flavour 
symmetry implementation, it is possible to reproduce a realistic sequential symmetry breaking 
chain, with the usual SM or MSSM Higgs fields at the electroweak scale, the introduction of a 
flavour symmetry puts strong constraints. Indeed we need additional scalars which transform 
under the gauge group and the effect is to sandwich the energy scales of the different symmetry 
breakings, lowering to 10 14 GeV the energy scale of the (almost) unification. 



3. Conclusions 

We have illustrated two models based on the flavour symmetry S4 x Z4 x U(1)fn where the 
BM mixing is realised at the LO in a natural way. The hierarchy of charged lepton masses is 
obtained as a combined effect of the U{1)fn and of S4 x Z4 symmetry breaking. 

Since exact BM mixing implies a value of tan #12 which is excluded by the data, large 
corrections are needed. The dominant corrections to the BM mixing arise at the NLO. The 
shifts of the quantities sin 2 9\2 and sin #13 from the BM values are linear in the parameter v' , 
which is expected to be of the same order as v, but not necessarily too close, as v and v' 
are determined by two different sets of minimisation equations. From the experimental value 
tan 2 9\2 = 0.45 =t 0.04, which is sizably different than the BM value tan 2 #12 = 1, we need 
v' ~ 0(A). As in most models where the BM mixing is only corrected by the effect of charged 



lepton diagonalisation, one also expects #13 ~ 0(A). A value of #13 near the present bound 
would be a strong indication in favour of this mechanism and a hint that the closeness of the 
measured values of the mixing angles to the TB values may be purely an accident. In addition, 
a very important feature of our models is that the shift of sin 2 #23 from the maximal mixing 
value of 1/2 vanishes at the NLO and is expected to be of 0(A 2 ) at most. In our 5*4 models, 
this property is obtained by only allowing the breaking of S4 in the neutrino sector via flavons 
transforming as 1 and 3 (in particular with no doublets). 

The quark sector is discussed only in the second model and realistic fermion mass hierarchies 
and a correct CKM matrix are found. This model is constructed in a supersymmetric Pati-Salam 
model, where a combined study of the flavour and Higgs sectors is performed. 

Studies on flavour changing neutral currents as in [241 [25 , 26J and on the Higgs phenomenology 
as in [271 [281 [29] will follow. 

Acknowledgments 

We thank the organizers of DISCRETE '10 - Symposium on Prospects in the Physics of Discrete 
Symmetries for giving the opportunity to present my talk and for the kind hospitality in Rome. 

References 

[1] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo, and A. M. Rotunno, in the proceedings of 4th International 
Workshop on Neutrino Oscillations in Venice: Ten Years after the Neutrino Oscillations, Venice, Italy, 
15-18 Apr 2008, pp. 21-28 [arXiv: 0809.2936]. 

[2] T. Schwetz, M. A. Tortola, and J. W. F. Valle, New J. Phys. 10 (2008) 113011 [arXiv: 0808.2016]. 

[3] P. F. Harrison, D. H. Perkins, and W. G. Scott, Phys. Lett. B530 (2002) 167 [arXiv: hep-ph/0202074]. 

[4] Z.-z. Xing, Phys. Lett. B533 (2002) 85-93 [arXiv: hep-ph/0204049] . 

[5] E. Ma and G. Rajasekaran, Phys. Rev. D64 (2001) 113012 [arXiv: hep-ph/0106291]. 

[6] K. S. Babu, E. Ma, and J. W. F. Valle, Phys. Lett. B552 (2003) 207-213 [arXiv: hep-ph/0206292]. 

[7] G. Altarelli and F. Feruglio, Nucl. Phys. B720 (2005) 64-88 [arXiv: hep-ph/0504165]. 

[8] G. Altarelli and F. Feruglio, Nucl. Phys. B741 (2006) 215-235 [arXiv: hep-ph/0512103]. 

[9] A. Y. Smirnov, Talk given at DISCRETETO - Symposium on Prospects in the Physics of Discrete Symmetries. 
[10] F. Feruglio, C. Hagedorn, Y. Lin, and L. Merlo, Nucl. Phys. B775 (2007) 120-142 [arXiv: hep-ph/0702194]. 
[11] F. Bazzocchi, L. Merlo, and S. Morisi, Nucl. Phys. B816 (2009) 204-226 [arXiv: 0901.2086]. 
[12] F. Bazzocchi, L. Merlo, and S. Morisi, Phys. Rev. D80 (2009) 053003 [arXiv: 0902.2849]. 
[13] I. de Medeiros Varzielas and G. G. Ross, Nucl. Phys. B733 (2006) 31-47 [arXiv: hep-ph/0507176]. 
[14] I. de Medeiros Varzielas, S. F. King, and G. G. Ross, Phys. Lett. B644 (2007) 153-157 [arXiv: hep- 
ph/0512313]. 

[15] Y. Lin, Nucl. Phys. B813 (2009) 91-105 [arXiv: 0804.2867]. 

[16] A. Y. Smirnov, Phys. Scripta T121 (2005) 57-64 [arXiv: hep-ph/0412391]. 

[17] M. Raidal, Phys. Rev. Lett. 93 (2004) 161801 [arXiv: hep-ph/0404046]. 

[18] H. Minakata and A. Y. Smirnov, Phys. Rev. D70 (2004) 073009 [arXiv: hep-ph/0405088]. 

[19] F. Vissani, arXiv: hep-ph/9708483. 

[20] G. Altarelli, F. Feruglio, and L. Merlo, JHEP 05 (2009) 020 [arXiv: 0903.1940]. 

[21] R. de Adelhart Toorop, F. Bazzocchi, and L. Merlo, JHEP 08 (2010) 001 [arXiv: 1003.4502]. 

[22] E. Bertuzzo, P. Di Bari, F. Feruglio, and E. Nardi, JHEP 11 (2009) 036 [arXiv: 0908.0161]. 

[23] D. Aristizabal Sierra, F. Bazzocchi, I. de Medeiros Varzielas, L. Merlo, and S. Morisi, Nucl. Phys. B827 

(2010) 34-58 [arXiv: 0908.0907]. 
[24] F. Feruglio, C. Hagedorn, Y. Lin, and L. Merlo, Nucl. Phys. B809 (2009) 218-243 [arXiv: 0807.3160]. 
[25] F. Feruglio, C. Hagedorn, Y. Lin, and L. Merlo, Nucl. Phys. B832 (2009) 251-288 [arXiv: 0911.3874]. 
[26] F. Feruglio, C. Hagedorn, and L. Merlo, JHEP 03 (2010) 084 [arXiv: 0910.4058]. 
[27] R. de Adelhart Toorop, F. Bazzocchi, L. Merlo, and A. Paris, arxiv: 1012.1791. 
[28] R. de Adelhart Toorop, F. Bazzocchi, L. Merlo, and A. Paris, arxiv: 1012.2091. 

[29] R. de Adelhart Toorop, Talk given at DISCRETETO - Symposium on Prospects in the Physics of Discrete 
Symmetries. 



